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MODAL DENSITY FUNCTION AND NUMBER OF 
PROPAGATING MODES IN DUCTS 
by E. J. Rice 
Lewis Research Center 

ABSTRACT 

Often raised questions in duct sound propagation studies involve the 
total number of propagating modes, the number of propagating radial 
modes for a particular spinning lobe number, and the number of modes 
possible between two given values of cutoff ratio or eigenvalue. These 
questions can be answered approximately by using the modal distribution 
function which is the integral of the modal density function for ducts in a 
manner similar to that previously published for architectural, acoustics. 
The modal density functions are derived for rectangular and circular 
ducts with a uniform steady flow. Results from this continuous theory 
are compared to the actual (discrete) modal distributions. 

INTRODUCTION 

The question of the number of propagating modes within a small 
range of mode cut-off ratio has been raised in references 1 and 2. The 
population density of modes has been shown to be greatest near cut -off 
and least for the well propagating modes. It has been shown in references 
1 and 2 that modes of nearly the same cut-off ratio behave nearly the same 
in a sound absorbing duct as well as in the way they propagate to the far 
field outside of the duct. Thus it has been proposed that, rather than 
handling all of the propagating modes individually, they can be grouped 
into several cut-off ratio ranges. It is important to know the modal den- 
sity function to estimate acoustic power distribution in the extension of the 
above mentioned studies. 
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A natural companion of the modal density function is the number of 
propagating modes possible in a duct. In the more limited case of rotor- 
stator interaction noise where the mode lobe numbers (m) are restricted 
(ref. 3) it is of interest to know the number of propagating radial modes 
for a particular mode lobe number. It is of course possible to enumerate 
the number of propagating modes by counting them with the use of eigen- 
value tables as done in reference 4, but it is much more convenient to 
have approximate equations for this purpose particularly when large num- 
bers of modes are involved. 

Equations for the number and modal density of the propagating modes 
in circular and rectangular ducts are developed in this paper. The basis 
for the equations lies in the very analogous problem of normal modes in 
rooms used in architectural acoustics (refs. 5, 6 are examples). 


SYMBOLS 

A„„ area of the axial modes in the j,m vector plane, see eq. (27) 
Aq B area of the oblique modes in the j,m vector plane, see eq. (29) 
c speed of sound, m/sec 

modal density function, see eqs. (40) and (42) 

D circular duet diameter, m 

F function of 77, £, and M, see eq. (26) 

Fq limiting value of F at a particular value of cut-off ratio 

f frequency, Hz 

H rectangular duct height, m 

i V-i 

j integer index, see eq. (25) 

J Bessel function of first kind of order m 

m 
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N total number of modes 

^ax num ^ er modes including 3=0 and m=0 modes only 

Nq number of modes with m=0 

N number of modes with a particular value of m 
N* number of modes, sum of N Q and N + 

N + number of spinning modes in one transverse direction 
M uniform steady flow Mach number 

m spinning mode lobe number 

maximum value of m for which first radial mode will propagate 

n 

P acoustic pressure, N/m 

R ratio of large to small dimensions in a rectangular duct 
r radial coordinate, m 

r^ hub radius in an annulus, m 

r w outer wall radius, m 

t time, sec 

W rectangular duct width, m 

x axial coordinate, m 

y transverse coordinate between rectangular duct walls, m 
z transverse coordinate, m 

or eigenvalue 

o!q limiting value of the eigenvalue for which propagation occurs 
/3 eigenvalue 

6 annulus hub -tip ratio, r^/r^ 

i] frequency parameter, fD/c or fH/c 

0 angular coordinate, radians 
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£ mode cut-off ratio, eq. (5) for cylindrical ducts, eq. (24) for 
thin rectangular ducts 

t propagation coefficient, see eqs. (1), (2), (20), (21), (35) 

Tq function of m and see eq, (8) 
u> circular frequency, rad/sec 


NUMBER 1 . OF PROPAGATING MODES 

Expressions for the number of propagating modes will be derived in 
this section. Both circular and rectangular ducts will be considered. 
These equations will be used in the next section to derive the modal den- 
sity functions. 


Circular Ducts 


The expression for the acoustic pressure of a spinning mode in a 
circular duct is given by 


where 


i wt-ini 0-i—Tx 


P = J 
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T = 


-M 4’ / / 1 - (1 - M 2 )(— j 
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1 - M 


and 
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c 


( 3 ) 


and M is the Mach number of the uniform steady flow in the duct. The 
condition that the radial particle velocity must vanish at the outer hard 
wall determines the eigenvalues from 


J>)=0 (4) 

The problem of modal enumeration amounts to counting the number of solu- 
tions to equation (4) such that a is less than some The quantity 
will be interpreted here as the limiting value of the eigenvalue for which 
sound propagation can occur. This idea of propagation involves the cut-off 
ratio defined as 

i- - -"L- (a 

ayi -M 2 

When equation (5) is used in equation (2), it is seen that at £ = 1 the radical 
in equation (2) vanishes and this will be termed cut-off. For £ < 1, the 
propagation coefficient t will be complex and damping occurs in the pres- 
sure (eq. (1)). Whether a mode is propagating or not depends upon the mode 
eigenvalue (a) as well as the duct diameter (D) ? sound frequency (f) and the 
uniform steady flow Mach number (M). Thus, for the purposes of this 
development, 


“0 - - 

V 1 -M^ 

The lobe number m^ is the largest value of m for which the first root of 
equation (4) is less than a q. 

The properties of the Bessel function eigenvalues will not be considered 
in detail here. For more detail, reference 5 is suggested and the following 
eigenvalue properties were obtained from this reference. The properties of 
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the eigenvalues are such that for a given lobe number m, the number of 
modes with eigenvalue a less than can be approximated by , 


a Q 1 

= — (sin T 0 ' T 0 cos T 0 ) + 4 
7T ^ 


CO 


where 


COS Tq = 


m 



( 8 ) 


Equation (7) is very good for large values of o>q and is fairly adequate 
even down to a q of ten. The total number of modes with a less than 
ckq is obtained by summing over equation (7) from m = 0 to uiq (see 
ref. 5). The result is, 


Qtr 


N*(a ft )- ii + 


Olr 


JT; 


0) 


(Note; there is an error in reference (5), eq. (15) in the last term). Equa- 
tion (9) contains the axisymmetric modes (m = 0) and the spinning modes in 
one transverse direction only. The number of axisymmetric modes is 
(eq. (7), m = 0), 


N 0 (c* 0 ) 


a. 


7T 


( 10 ) 


The number of spinning modes in one transverse direction is thus, 


N + (o- 0 ) = N* (a Q ) -_N 0 (a 0 ) 
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Finally the total number of modes in the positive axial direction is 

N(t» 0 ) = 2N + (a 0 ) + N 0 (a 0 ) = ^ + — (12) 

4 2 

The number of propagating modes is now related to the duct size, 
steady flow Mach number and sound frequency by using equation (6) to 
yield, 


and 



(13) 


(14) 


n - Z ' f * n — (is) 

4(1 ~M 2 ) 2Vl -M 2 ) 

The number of propagating modes is seen to depend only upon the fre- 
quency parameter 77 and the steady flow Mach number M„ The number of 
modes (spinning only, in. one direction) have been counted in reference 4 for 
M = 0, tj = 30. 11 with a reported result of N + = 1126. Equation (14) yields 
N + = 1127 showing excellent agreement for large values of 77 . For small 
values of 77 and for M = 0, figure 1 shows a comparison between an exact 
mode count and N from equation (9) (using eq. (6) with M = 0). The 
agreement is seen to be quite good. Note that the equations do not include 
the plane wave mode since when 77 == 0 then N = 0). If this mode were also 
included (add unity to equations) the agreement between calculated and actual 
mode counts would be improved at low values of 77 . 
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A more general form of equations (13) to (15) can be developed for an 
arbitrary value of the cut-off ratio rather than just at cut-off (£ = 1). This 
is done by using equation (5) rather than (6) for a or a. q. The results are, 


N 0 (5 )=— JL_ (16) 

i'fT. M 2 

2 2 

N + ({ ) = Z-2 + Y (17) 

8{ 2 (1-M 2 ) 2JA/i -M 2 

and 

N(£) = 1 Ljt + M .. - (18) 

4£ 2 (1-M 2 ) 2£ \fl~ M 2 

These give the number of propagating modes with cut-off ratios from infinity 
to the selected value of Equations (13) to (15) are thus special cases of 
equations (16) to (18) with £ = 1. 

The number of propagating modes with a particular lobe number (m), 
frequency parameter (77), and Mach number (M) may be of interest in a cir- 
cular duct. This number can be calculated directly from equation (7) using 
equations (8) and (6). An approximation to equation (7) can also be derived 
by using a series expansion around m = 0. The resulting approximation is 

N m ( m »’). M ) w + m V -. (19) 

4 2 2 „ 2 „ 

Sample calculations using equations (7) and (19) are compared to exact 
counts in table I for 77 = 15 and M = 0. 
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TABLE I. - COMPARISON OF EXACT AND CALCULATED N m 


m 

counted 

N m 

(eq. (7)) 

N m 

(eq. (19)) 

0 

15 

15. 25 

15. 25 

5 

13 

12. 83 

12. 83 

10 

11 

10. 59 

10. 59 

20 

7 

o. 62 

6. 60 

30 

3 

3.41 

3. 29 

35 

2 

2. 12 

1. 89 

40 

1 

1. 09 

0. 65 


The agreement using either equation is seen to be quite good. In actual 
use the calculated N would probably be truncated or rounded off to 
obtain an integer. 


Thin Rectangular Duct 

The thin rectangular duct is of interest since it approximates an un- 
wrapped thin annular duct (ref. 7) often used in noise suppressors. The 
boundary condition in the unterminated direction is that the pressure must 
repeat over a distance of one annulus circumference. The acoustic pres- 
sure solution is, 

• ,,,^- imz Jcutx 

P = cos e Tw (20) 

for symmetric modes, with sin (2o!y/H) replacing the cosine for the asym- 
metric modes. The axial propagation coefficient is 
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where now 



c 


and the hub -tip ratio 6 is, 


( 21 ) 


( 22 ) 



(23) 


Note that as in the previous section traveling waves in the positive z direc- 
tion only are included. 

The cut-off ratio is defined as, 



The hardwall condition at y = ±H/2 requires that 

a=VL, j = 0,l,2,3,... (25) 

2 

which includes the eigenvalues for both the sine and cosine functions. Equa- 
tion (24) can be rearranged (and using equation (25)) to yield, 



.1 •• M- 
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The enumeration of the modes reduces to finding the number of j and m 
values such that the value of F calculated from equation (26) is less than 
or equal to a limiting value Fq, This is the same as counting the modes 
with cut-off ratios greater or equal to a particular value of £„ Due to the 
even spacing of the eigenvalues it is easiest to use the vector plane approach 
of reference 6. Equation (26) may be viewed as a vector equation in the j 
and m plane and is pictured as such in figure 2. The legation of each mode 
is shown by a dot and has associated with it a unit area of (1 - 6)/jr. The 
term axial modes refer to those along both the j = 0 and m = 0 axes while 
the oblique modes are all of the remaining modes. The idea is to calculate 
the area associated with each type of mode and divide this by the unit area to 
obtain the number of modes. The area associated with the axial modes is 



1 + 


(1 - 5 ) 


IT J 


(27) 


that part being associated with m = 0 being pictured in figure 2. The num- 
ber of axial modes is thus , 


N ax = F 0 


1 + 


7T 


(1 " 6)J 


(28) 


The area for the oblique modes is 


rrF? F 


l OB 


o.io ! . (i - a) 

2 L n - 


(29) 


where half of the area associated with the axial modes has been removed from 
the quarter circle since this area has already been accounted for with the axial 
modes. The final result is 
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N*(|) 


2 2 
7T 7? 


V 


| 2 (1 ~ M 2 )(1 - 6) 


7T 

(1 - 6)_ 


(30) 


which represents the m = 0 modes and the transverse modes in one direc- 
tion only. The number of modes associated with m = 0 (use eq. (26)) is, 


N o«> - -7= 

1 Vi - ^ 


(31) 


with spinning modes in one direction, 


N + (?) = 


iV 


V 


| 2 (1 -M 2 )(l - 6) £ Vl - m2 

and with, the total number of modes being 


trf- 


- 

1 

6) J 


(32) 


N((j) = 


2m] 


7TV 


£ Vl~- M 2 (1 - 6) 


vT 


M" 


.+ 1 


J 


(33) 


which represents the m = 0 modes and the modes spinning in (lie plus and 
minus transverse directions. The number of propagating modes of each 
type can now be determined by letting the cut-off ratio £ be unity. Figure 3 
shows a comparison of N (eq. (30)) with an actual count of modes using 
equation (25). No attempt has been made to reproduce the discrete steps of 
the actual count (as in. fig. I). The agreement is seen to he quite good and 
thus the continuous equations can be used to approximate the actual discrete 
number of modes. Included for comparison is a mode count (using tables of 
ref. 8) in an annulus with hub-tip ratio (8) of 2/3, The actual annulus is seen 
to have somewhat fewer modes than its approximating rectangle but the agree- 
ment should improve for higher 8. 

The number of modes propagating with a particular lobe number is quite 
simple to obtain, with a thin rectangular duct. Solve for j in equation (26) 
(with £ = 1) and add unity (to account for the j = 0 mode). The result is, 
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N m (m } ??,M) = 


/ 4t? 2 _ 

~m(l - 6)' 

/ (1 - M 2 ) 

7T 


+ 1 


(34) 


Note that equations (30), (31), and (33) do not contain the plane wave mode. 
If the plane wave is desired in the count, unity should be added to these 
equations. 


Rectangular Duct 

Consider a rectangular duct with smaller dimension H and larger di- 
mension W. The pressure for the symmetric modes is given by 


P = cos 




• . -CO 

io)t-i-rx 


H 


W 


(35) 


with sine functions used for the asymmetrical modes. A procedure very 
similar to that in the previous section will yield, 

N($> = — — . + JlQ LtlL (36) 

£ 2 (1 - M 2 ) I t\fl - M 2 

with 


and 



c 



.. // 

// 


R=W 

H 


> 1 


(37) 


(38) 
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Note the frequency parameter is based upon the smaller dimension here. 

In this case no arguments involving plus or minus direction transverse 
waves need be made since the cosine or sine function encompass both wave 
directions. Again with £ = 1 in equation (34) the total number of prop- 
agating modes can be estimated. 

Note that all of the mode counting equations previously presented con- 
sider only the modes propagating in the plux-x direction. If reflections 
off of an end termination are considered, the number of modes would be 
doubled. 


MODAL DENSITY FUNCTION 

The modal density function for all modes in a circular duct will be 
developed here. Similar procedures can be used for the axisymmetric or 
the spinning modes or for the other duct shapes if desired. 

Differentiation of equation (18) yields 


dN(j) . jV 

2J 3 (1-M 2 ) 

where the sign has been changed such that integration from £ to °° ■will 
recover the original mode number equation. When equation (39) is normal- 
ized by the total number of propagating modes (eq. (15)) the modal density 
function is obtained as 


1 + 






TTH] 


(39) 



2 

1 dN(£) 

i + Vi-m 2 4 

7T7] 

N d* 

1 + Wl -M 2 

nr) 


( 40 ) 
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The second term in the denominator is usually small for aircraft inlet appli- 
cations since usually 10 < 77 < 30. The second term in the numerator can 
dominate only if 

| . (41) 

Yl -M 2 

which is a very large cut-off ratio (far from cut-off (£ = 1) which is of most 
interest for acoustic power considerations). The modal density function can 
thus be well approximated by 

— («) 

The same equation can be used for rectangular ducts provided that 

g < - -JffL— (43) 

yi -m 2 

for thin rectangular ducts and 

£ < 2 mjR ( 44 ) 

f ]/ 1 - M 2 (R + 1) 

for rectangular ducts. Some care should be exercised in using equation (42) 
for rectangular ducts. In some applications conditions (43) or (44) may not 
be realized. In aircraft applications where annular ducts are used (rec- 
tangular intended to simulate annulus) the frequency parameter (?j) is usually 
around one or. two representing a limiting cut-off value beyond the range of 
interest. 

The modal density is seen to be greatest near cut-off (£ = 1) and dim- 
inishes rapidly as cut-off ratio is increased. This quantifies the observa- 
tion of this fact which was stated in reference (1). It is also interesting to 
note that the modal density function (in its approximate form) depends only 
upon the cut-off ratio. 
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CONCLUDING REMARKS 

An approximate modal density function has been derived which is valid 
for both circular and rectangular ducts. This function represents the frac- 
tion of the total modal population which is within a small cut-off ratio range. 
The modal density function was found to depend only upon cut-off ratio. 

This coincides with the approximate dependence upon cut-off ratio (refs. 1 
and 2) of optimum wall impedance, maximum possible sound attenuation, 
and far field sound radiation. The density function is required to extend an 
approximate liner design technique based upon the unifying influence of cut- 
off ratio (ref. 2). 

Approximate equations were also presented to estimate the number of 
propagating modes in circular and rectangular duets. All of the above are 
approximate in the sense that continuous equations are presented to repre- 
sent an inherently discrete phenomenon. The approximation is good pro- 
vided a sufficiently large number of modes can exist. In cases where only 
a few modes exist due to low frequency, small ducts, or restricted noise 
sources (such as some cases of rotor -stator interaction) then the sound 
propagation should be handled on an individual mode basis. 
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FROM EXACT MODE COUNT 



Figure 1. - Comparison of mode count with approximating equation for 
low values of frequency parameter. 


AXIAL MODE AREA FOR m>0 



Figure 2. - Vector plane area representation used 
for mode summation in the thin rectangular 
duct. 
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Figure 3. - Comparison ol mode counts with approximating 
equation, rectangular duct simulation of annulus. 
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